Critical Collapse of a Complex Scalar Field with Angular Momentum 
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We report a new critical solution found at the threshold of axisymmetric gravitational collapse 
of a complex scalar field with angular momentum. To carry angular momentum the scalar field 
cannot be axisymmetric; however, its azimuthal dependence is defined so that the resulting stress 
energy tensor and spacetime metric are axisymmetric. The critical solution found is non-spherical, 
discretely self-similar with an echoing exponent A = 0.42 ± 4%, and exhibits a scaling exponent 
7 = 0.11 ± 10% in near critical collapse. Our simulations suggest that the solution is universal 
(within the imposed symmetry class), modulo a family-dependent constant phase in the complex 
plane. 



The main purpose of this work is to study the effect 
of angular momentum in axisymmetric critical collapse 
of massless scalar fields. Critical collapse refers to the 
threshold of black hole formation, where interesting ef- 
fects known as critical phenomenal have been observed 
in the gravitational collapse of a wide variety of types of 
matter, as well as vacuum gravitational waves (see 
for recent reviews). For spherically symmetric massless 
scalar collapse, this behavior includes universality, scale 
invariance, and power law scaling of length scales that 
arise near criticality. In super-critical collapse, the char- 
acteristic length is the mass, M, of black holes that form. 
In the case of rotating collapse, since angular momentum 
has dimension length , one might naively expect the an- 
gular momentum, J, of the black holes formed to scale 
as J oc M^. A more refined analysis carried out using 
perturbation theory |3] suggests that J (x Af ^(^"^"^[^'l) , 
where Re[A'] is the real part of the exponent A' of the 
dominant perturbative mode that carries angular mo- 
mentum. In i. Re [A'] was found to be roughly —0.017, 
implying an approximate scaling J oc M^ °^. Thus, the 
Kerr parameter a = J/M'^ is expected to scale to zero 
(albeit slowly) as the black hole threshold is approached. 
Incidentally, this picture could be altered significantly if 
the numerical evidence of ^ proves correct, and there is 
a non-spherical unstable mode of the critical solution. 



Sec. m (see for an alternative approach). One draw- 
back to this method is that imposing such an ansatz for 
the complex field forces a non-spherical energy distri- 
bution. This means that the class of solutions we can 
study occupy a region of phase space distinct from that 
of spherical spacetimes, and so we cannot explore the 
role of angular momentum as a perturbation in spheri- 
cal critical collapse 0, 0- On a positive note, the fact 
that we do find a new (axisymmetric) critical solution is 
interesting aside from questions of angular momentum, 
because it suggests that phase space has a more intricate 
structure than one might have naively imagined, proba- 
bly containing an infinite set of distinct intermediate at- 
tractors characterized by their behavior near the center 
of symmetry (the results of are also in accord with this 
conjecture). Regarding the question of how net angular 
momentum affects threshold behavior in this model — it 
appears to be irrelevant, with the angular momentum 
of black holes formed in super-critical collapse decaying 
significantly faster than J cx Af^, though our simulations 
are not accurate enough to calculate exactly how fast. 
Below we briefly describe the physical system and nu- 
merical code we use (for more details see |E IE IS ITof'l- 
and present our results in Sec. ^1 

I. PHYSICAL SYSTEM 



In general, numerical exploration of angular momen- 
tum in the collapse of a single real scalar field would 
require a 3D code, for axisymmetric distributions of such 
a matter source cannot carry angular momentum. Con- 
structing a general relativistic 3D simulation capable of 
resolving the range of length scales that unfold in scalar 
field critical collapse is a daunting project, and may re- 
quire computational capacity not currently available. A 
"cheaper" alternative is to consider a set of distinct scalar 
fields, each with azimuthal dependence and hence an- 
gular momentum, and then add the different fields co- 
herently such that the net stress energy tensor is ax- 
isymmetric. A natural way to achieve such a coherent 
sum is via a single complex field, as will be explained in 



We consider the Einstein field equations 
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(1) 



where g^i/ is the spacetime metric, i?^^ is the Ricci ten- 
sor, R = Rf^fj^ is the Ricci scalar, and we use geometric 
units with Newton's constant G and the speed of light c 
set to 1. We use a massless, minimally-coupled, complex 
scalar field '5 (with complex conjugate 4") as the matter 
source. 4* satisfies a wave equation 4':^'' = 0, and has a 
stress energy tensor T^^, given by 



(2) 
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We solve Q and the wave equation (hereafter the field 
equations) in axisymmetry, using coordinates [t,p,z,(j)], 
where (j) is adapted to the azimuthal symmetry, t is time- 
like, and (p, z) reduce to standard cylindrical coordinates 
in the flatspace limit. The axial Killing vector is then 



(3) 



The existence of this Killing vector allows us to define 
the conserved angular momentum, J, of the spacetime 



J = - 



(4) 



where the integration is over the t = const, spacelike hy- 
persurface, S, h is the determinant of the intrinsic metric 
on E, and is the hypersurface normal vector. Using 
© and Q, O evaluates to 



J 



(5) 



Thus, for a configuration of the scalar field to have non- 
zero angular momentum, ^' must have some azimuthal 
dependence. We thus adopt the following ansatz 



^{p,z,t,(f>) = ^{p,z,t)e 



imcj) 



(6) 



where $(p, z,i) is complex, and m must be an integer 
for the scalar field to be regular. It is straightforward to 
check that this form of 4" gives a stress-energy tensor that 
is 0-independent, yet can yield net angular momentum. 
Note that the on-axis {p = 0) regularity condition for 
<I> depends upon the value of m; specifically we must 
have \imp^o^{p, z,t) = p'^f{z,t). For simplicity and 
specificity, we hereafter restrict attention to the case m = 
1. As stated above, we expect that additional, distinct 
critical solutions exist for m — 2,3, ■ ■ ■. 

To keep our discussion concise, we only state the metric 
and set of variables we use (all functions of (p, z,t)), and 
briefiy describe the solution procedure — more details can 
be found in IM la la Un. . The line element is 



ds^ = -a^dt^ + (7) 
i;^[{dp + pPdtf + (dz + P'-dtf + p2e2p*d02] _^ 

The lapse function, a, is fixed by maximal slicing, and (Uj) 
refiects the additional coordinate conditions we have im- 
posed: conformal fiatness of the two dimensional p — z 
subspace, and = 0. The Einstein equations are written 
in first-order-in-time form by introducing the following 
"conjugate" variables 

n = {-2K/ - K/') Ip, uj^ = e^p^s£,f'e'\ (8) 
where is the extrinsic curvature tensor and is 



the "twist" of the Killing vector. We separately evolve 
the real and imaginary components of the scalar field by 
defining real functions and 3>i via 

^ = p{^r + ^<^i) (9) 

and their dynamical conjugates H,. and 11^ by 

= Re[$^an"]/p, n, = Im[$,„77"]/p. (10) 

The factors of p appearing in the above definitions are 
included so that ^i, and 11^ satisfy Neumann con- 
ditions on-axis. Similarly, the variables corresponding to 
LJa and that are evolved in the code have appropriate 
powers of p factored out so that they satisfy Dirichlet 
conditions on-axis (see ^3 for the specific definitions). 

For the evolutions presented here we use the following 
initial data for the scalar field components 



$ru(p, z,0) = ^rljexp 



Iir\t{P-,Z,Q) = er\i<^r\i{P,Z,{)) , 



- p^ + Z^ - Rr\A /S 



(11) 



where Ar,Ai,Rr,Ri,Y,,€r and ti are parameters fixing the 
shape of the initial scalar field profiles. All other freely 
specifiable variables are set to zero at i = 0, while the 
constrained variables a,^ and {j3P , (3^ ,Wp\ arc obtained 
by solving the maximal slicing condition, Hamiltonian 
and momentum constraints respectively. 

We use a partially constrained finite-difference scheme 
with adaptive mesh refinement (AMR) to evolve the sys- 
tem of equations with time. In particular, the slicing con- 
dition and momentum constraints are used to fix a and 
{(]'', (3^} respectively, while the remainder of the variables 
are updated using their evolution equations. 

II. RESULTS 

We now present results from a preliminary study of 
the black hole threshold of the complex scalar field sys- 
tem introduced in the previous section. We focus on four 
sets of initial data, summarized in Table Family A is 
the 'canonical' example, consisting of identical pulses of 
^riPi z, 0) and <&i(p, z, 0) that are initially approximately 
outgoing {tr — 1) and ingoing (ej = —1) respectively. 
This choice for (6^,6^) in a sense maximizes the net an- 
gular momentum (jSJ , given the initial profiles for and 
^i. Conversely, family B is time-symmetric, and hence 
has zero net angular momentum. Families A and B can 
be written as <i>(p, 2,0) = A{p, z)e'^^° , where A[p,z) is a 
real amplitude function and 5q is a constant phase factor, 
equal to 7r/4 in both cases. Family C is thus identical to 
family A except for the initial phase. For family D 
and have distinct initial profiles and thus cannot be 
characterized by a constant phase. 

Based upon the collapse simulations we have per- 
formed for these four families of initial data, we can 
suggest the following about the threshold behavior for 
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Label 


P 


Rr 


Ri 




e,. 


e-i 


5o 


5* 


A 


Ar = Ai 


0.6 


0.6 


0.1 


1 


-1 


7r/4 


0.91 ±3% 


B 


A-r ^ Ai 


0.6 


0.6 


0.1 








7r/4 


7r/4 


C 


Arf ^ S-Zl^ 


0.6 


0.6 


0.1 


1 


-1 




0.39 ±3% 


D 


Ar = Ai 


0.65 


0.6 


0.1 


1 


-1 




1.34 ±3% 



TABLE I: Parameters (see pi^ l for the four families of initial 
data discussed here, where p denotes the parameter(s) we tune 
when searching for the black hole threshold. 5o and 5* are the 
phase of the initial data (if applicable) and estimated phase of 
the critical solution respectively. All of the simulations were 
run with the outer boundary (pmax, Zmax, —z^in) at 4 (except 
for the data of Fig. 2] where the outer boundary was at 16), 
and in each case p was tuned to within {p — p*)/p* « lO"^** 
of threshold. 




FIG. 1: A surface plot of the real part {'^rp) of the com- 
plex field after several echoes of a near-critical evolution. The 
origin is at the bottom-center of the figure, the p axis runs 
vertically through the middle, and the z axis runs horizon- 
tally. Only a single echo (roughly) at the origin corresponds 
to the self-similar part of the spacetime — the other "waves" 
visible were radiated during earlier echoes of the field. Note 
also that the solution is not spherically symmetric. 



this matter model. There is apparently a discretely 
self-similar critical solution that is universal to within 
a family-dependent phase. In other words, one can write 
the critical solution $* for the scalar field as $*(p, z,t) = 
A*{p, z,t)e^^ , where A*{p,z,t) is a universal real func- 
tion and S* is a family dependent constant (see Table 
U). Note that this phase dependence is a consequence 
of the U{1) symmetry of the Lagrangian of the complex 
field, and has been observed in charged scalar field crit- 
ical collapse^2- Also, note that any self-similar solu- 
tion is unique only up to a global rescaling of the form 
(f, ic*) (nt^Kx'') when written in suitable coordinates 
{t, a;*), with k a constant. Fig.^shows a snapshot of the 
real part of the scalar field ($rp) at late times in a near- 
critical collapse simulation. To estimate S* for a given 
family, and the echoing exponent A for the putative crit- 
ical solution, we examine the central value of the real 
and imaginary parts of the scalar field divided by proper 
radius pc (to factor out the leading order approach to 
of $ in a covariant manner) 

Pe^pV'V'^. (12) 

Fig. 121 shows plots of $,.rp//3c and ^ivp/pc vs. — Inr 
for the nearest-to-threshold solutions found, where r is 
central proper time (see Q for our definitions and compu- 
tations of quantities measured by central observers, but 
note that we define t such that the accumulation event 
of the critical solution corresponds to r = 0). We have 
multiplied the scalar field by r in these plots to cancel the 
artificial growth introduced by dividing by pc- Note that 
the equations of motion for $r and $i are identical; hence 
if <&r(p, 0) = ^'i(p, 2, 0) (as with family B), then the ini- 
tial phase, (5o = 7r/4, is preserved during evolution. The 
echoing exponent A is the period of the self-similar solu- 
tion in logarithmic proper time; from Fig. |21 (and similar 
data for family B) we estimate A = 0.42 ± 4%. To esti- 
mate the scaling exponent 7, we measure how the maxi- 
mum value attained by the Ricci scalar (on axis), in 
sub-critical evolutions depends upon the parameter-space 
distance from threshold, p* — p 0|. Representative re- 



r ^ rp/p^ r|'.Tp/p 




J \ \ L 



1 2 3 4 5 

— IriT 

FIG. 2: The real (p$r) and imaginary (p$r-) components of 
the central value of $ @ multiplied by proper time r and 
divided by proper radius pc 112II . versus — Inr for the near- 
critical collapse solutions of families A, C and D (the phase 
information for family B is trivial as <!>,. = <I>i then, and so 
for brevity we do not show it). The family D solution shown 
here is super-critical, and the simulation is stopped soon after 
an apparent horizon is detected. 



suits are shown in Fig. 13 Combining such data from all 
the families, we estimate 7 — 0.11±10%. For a discretely 
self-similar solution, one expects the linear relationship 
assumed in Fig . El to be modulated by an oscillation of 
period 2 A Il3l |: we have not run a sufficient number of 
simulations to adequately resolve such an oscillation. 

The uncertainties quoted above for 7, A and 5* (in Ta- 
ble P) were estimated from convergence calculations from 
simulations using 3 different values of the maximum trun- 
cation error threshold that controls the AMR algorithm, 
but do not account for possible systematic errors (see a 
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relationship closer to Jah « Mah^; however, we are still 
far from threshold there, and furthermore are entering 
the regime where the angular momentum calculation is 
dominated by numerical errors. Thus, perhaps the only 
quantitative statement we can make regarding angular 
momentum scaling for this system is that J scales to 
zero significantly faster that J cx M^. 



FIG. 3: Rm, the maximum value of the Ricci scalar on axis 
(p = 0) attained during sub-critical evolution, versus distance 
\p\ = p* ~ p from threshold, for family A. Each point repre- 
sents a single simulation. The line is a least-squares fit to the 
data; Rm has dimension length'^, and hence the slope of the 
fit is expected to be ~2"f. For this case we infer 7 « O.f 1. 

discussion of related issues in |^). 

Finally, as mentioned in the introduction, for the near- 
critical solutions described here, net angular momentum 
seems to be completely irrelevant. To within the accuracy 
of our simulations, we cannot differentiate between the 
late time, self-similar regions of the spacetimes obtained 
from families A or B, and in the latter case, there is 
no angular momentum. In fact, any angular momentum 
present is radiated away so rapidly during self-similar col- 
lapse that we cannot accurately calculate the correspond- 
ing scaling exponent (i.e. the remaining angular momen- 
tum is zero to within truncation error) . Fig. ^ shows a 
plot of the mass estimate Mah versus angular momen- 
tum Jah, on a logarithmic scale, of black holes formed 
in super-critical collapse. AIah and Jah are calculated 
from the area and angular momentum of the apparent 
horizon respectively (using the dynamical horizon frame- 
work ^3), and are computed at the time the apparent 
horizon is first detected. For large black holes (i.e. those 
with Mah of order the ADM mass). Fig. 0] suggests that 
Jah oc Mah^- However, this region of parameter space is 
"maximally" far from threshold, in that these are almost 
the largest black holes that we can form from initial data 
not already containing an apparent horizon. For some- 
what smaller black holes. Fig. 0] shows a transition to a 
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FIG. 4: Estimated black hole mass (Mah) versus angular 
momentum (Jah) in super-critical collapse of family A initial 
data. Points represent individual simulations, while the two 
lines are separate linear regression fits to the set of points to 
the left and right of the "knee" in the curve at In(MAH) ~ 
—3.2, with slopes ~ 6.0 and ~ 2.2 respectively. In Inp, the 
horizontal scale ranges from —22 to —2 (compare to Fig.j^J. 



[1] M.W. Choptuik, Phys. Rev. Lett. 70, 9 (1993) 

[2] C. Gundlach, Phys. Rept, 376, 339 (2003); C. Gundlach, 

Living Rev.Rel. 2 1999-4; A. Wang, Braz.J.Phys. 31, 188 

(2001) 

[3] D. Garfinkle, C. Gundlach, and J. M. Martin-Garcia 

Phys. Rev. D59, 104012 (1999) 
[4] J. M. Marti'n-Garci'a and C. Gundlach, Phys. Rev. D59, 

064031 (1999) 

[5] M.W. Choptuik, E.W. Hirschmann, S.L. Liebling and F. 

Pretorius, Phys. Rev., D 68 044007, (2003). 
[6] M.W. Choptuik, I. Olabarrieta, W.G. Unruh and J.F. 

Ventrella, in preparation 
[7] E.W. Hirschmann and D.M. Eardley, Phys. Rev. D51 

4198, (1995) 

[8] M.W. Choptuik, E.W. Hirschmann, S.L. Liebling and F. 



Pretorius, Class. Quant. Grav 20, 1857 (2003) 
[9] F. Pretorius, "Adaptive Mesh Refinement for Coupled 

Elliptic- Hyperbolic Systems", in preparation. 
[10] M.W. Choptuik, E.W. Hirschmann, S.L. Liebling and F. 

Pretorius, in preparation 
[11] C. Gundlach and J.M. Martin-Garcia, Phys. Rev. D 54, 

7353 (1996); S. Hod and T. Piran, Phys. Rev. D 55, 3485 

(1996) 

[12] D. Garfinkle and G.C. Duncan, Phys. Rev. D58, 064024 
(1998) 

[13] C. Gundlach, Phys. Rev. D55, 695 (1997); S. Hod and T. 

Piran, Phys. Rev. D 55, 440 (1997) 
[14] A. Ashtekar and B. Krishnan, Phys. Rev. Lett. 89 261101, 

(2002). 



